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The mechanical properties of light have found widespread use in the manipulation of gas-phase 
atoms and ions, helping create new states of matter and realize complex quantum interactions. 
The held of cavity-optomechanics strives to scale this interaction to much larger, even human-sized 
mechanical objects. Going beyond the canonical Fabry-Perot cavity with a movable mirror, here 
we explore a new paradigm in which multiple cavity-optomechanical elements are wired together to 
form optomechanical circuits. Using a pair of optomechanical cavities coupled together via a phonon 
waveguide we demonstrate a tunable delay and hlter for microwave-over-optical signal processing. In 
addition, we realize a tight-binding form of mechanical coupling between distant optomechanical cav¬ 
ities, leading to direct phonon exchange without dissipation in the waveguide. These measurements 
indicate the feasibility of phonon-routing based information processing in optomechanical crystal 
circuitry, and further, to the possibility of realizing topological phases of photons and phonons in 
optomechanical cavity lattices. 


I. INTRODUCTION 

Microscopic optical and microwave cavities, with their 
wavelength or even sub-wavelength mode size, have re¬ 
cently been used to greatly enhance the radiation pres¬ 
sure interaction between electromagnetic waves trapped 
in the cavity and mechanical vibrations of the cavity 
walls [T]. Technical advances in our ability to fabricate 
structures of ever smaller size and higher quality has led 
to the demonstration of a number of new phenomena, in¬ 
cluding laser cooling of mechanical resonators into their 
quantum ground state 0131 , optomechanically-induced 
electromagnetic transparency im, and microwave sig¬ 
nal processing 0 0. To date these experiments have 
largely involved single cavity elements. A new paradigm, 
enabled by chip-scale cavity-optomechanics, involves the 
integration of multiple optomechanical cavities together 
to realize extended photon-phonon excitations for classi¬ 
cal and quantum information processing applications 0 - 
m- Among the different types of optomechanical sys¬ 
tems currently being studied, the device architecture 
based on thin-film optomechanical crystals is 

a particularly promising integration platform given the 
large attainable radiation-pressure coupling and the abil¬ 
ity to create photonic and phononic band-gap waveg¬ 
uides. 

In an integrated optomechanical network composed of 
cavities and waveguides, propagating phonons can be 
routed and stored among different nodes as controlled 
by driving laser fields, providing various information pro¬ 
cessing functions - such as buffering, delay, and filtering - 
in a chip-scale footprint due to the slow velocity of acous¬ 
tic waves [ME]- As an initial step towards realizing 
more complex multi-element cavity-optomechanical sys¬ 
tems, here we explore the optical excitation and routing 
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of GHz phonons in a small optomechanical circuit con¬ 
sisting of two optomechanical crystal cavities connected 
by a dispersion-engineered phonon waveguide. Pulsed 
and continuous-wave measurements are first used to char¬ 
acterize the coupling efficiency and phonon propagation 
properties of the system. Utilizing separate optical driv¬ 
ing fields for each optomechanical cavity, we use this sys¬ 
tem to demonstrate an optically-tunable microwave de¬ 
lay line in which microwave-over-optical signals sent in 
one cavity port are efficiently converted into outgoing 
microwave-over-optical signals on the other cavity port. 


Owing to the continuum of modes supported by the 
phonon waveguide, a tight-binding type of interaction 
between distant cavities may also be realized. Such ideas 
have recently been explored utilizing photonic waveg¬ 
uides, both in the optical domain with photonic crystal 
cavities [20] and in the microwave domain with super¬ 
conducting qubits [21]. Here we use a phonon waveg¬ 
uide to strongly couple two distant optomechanical cav¬ 
ities, leading to an exchange-type phonon routing. The 
hybridized mechanical cavity modes, while spatially ex¬ 
tended, nevertheless have strong optomechanical cou¬ 
pling with both localized optical cavity modes which are 
not hybridized. Such a feature can be used to create 
mechanical transducers capable of bridging two sepa¬ 
rated physical systems [22], such as in quantum wave¬ 
length conversion between optical and superconducting 
microwave systems [231126] and coupling between solid 
state spins [T7ll29] . Scaling such strong coupling to arrays 
of cavity elements may also enable novel optomechani¬ 
cal metamaterials, where in combination with a spatially 
distributed phase pattern of the optical driving field m, 
topological phases for photons and phonons are expected 
to arise miEa. 
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FIG. 1. GHz dispersion-engineerable phonon waveguide, a, A schematic diagram of the single cavity-waveguide device. 
a and b represent the optical and mechanical cavity modes. The mechanical cavity mode couples to the phonon waveguide at 
rate 7 e and has an optically induced decay rate 70M. b, Scanning electron microscope (SEM) image of a silicon optomechanical 
cavity-waveguide device based on the scheme of a. c, Theoretical band structure of the phonon waveguide employed in this 
work. Green curve is the waveguide breathing mode and yellow curves are the bands of the phononic shielding. Inset on the 
right is the modal prohle of the waveguide breathing mode corresponding to wavevector /c = 0.8^. d, Optically transduced 
mechanical power spectral density (PSD) of a typical device showing a few hybridized cavity-waveguide resonances, e, Group 
velocity as a function of frequency across the breathing mode waveguide band. Green is the theoretical curve and red dots are 
the measured data. 


II. PHONON WAVEGUIDE 
CHARACTERIZATION 

While integrated phonon waveguides with frequency 
bands of kHz to a few hundred MHz have been pre¬ 
viously studied j33H35] . here we explore the properties 
of thin-film, dispersion engineered phonon waveguides 
in the GHz frequency range. A schematic diagram and 
SEM image of a cavity-waveguide device used to study 
such phonon waveguide properties are shown in Fig. 
and Fig. [^, respectively. The device is fabricated from 
the 220 nm silicon device layer of a silicon-on-insulator 
wafer (see App. 0 , and consists of several parts: an op¬ 
tomechanical cavity with co-localized optical and acous¬ 
tic resonances, an optical coupling waveguide, and a long 
phonon waveguide section which is end-fire coupled to the 
cavity. The optomechanical cavity is realized in a pho¬ 
tonic crystal nanobeam which supports an optical reso¬ 
nance at Ac ~ 1550 nm in the telecom wavelength band 
that is strongly coupled through the elasto-optic radi¬ 
ation pressure effect to a “breathing mode” mechanical 
resonance of the beam at a frequency of co’^/27r ~ 6 GHz. 
One end of the nanobeam is terminated by a two di¬ 
mensional phononic crystal mirror with a full phononic 
bandgap in the relevant frequency range m-, the other 
end is connected to the phonon waveguide consisting of 


a similarly patterned central beam region surrounded by 
phononic bandgap shielding on either side. The phonon 
waveguide is designed to have large extrinsic coupling 
7e to the nanobeam cavity mode, leading to hybridiza¬ 
tion of the localized cavity mode and the nearly-resonant 
phonon waveguide modes. The optical coupling waveg¬ 
uide is fabricated in the near-field of the nanobeam cavity 
(bottom beam of Fig. Ef)- allowing for evanescent cou¬ 
pling of laser light into and out of the cavity. A single 
optical fiber taper is used to couple light into the on- 
chip coupling waveguide, and a photonic crystal mirror 
is etched in to the end of the optical coupling waveguide 
so that light coupled into the nanobeam cavity can be 
recollected by the optical fiber taper as per Ref. [36] . 

Figure shows the band structure of the phonon 
waveguide with silicon slab thickness d = 220 nm and 
lattice constant a = uq = 480 nm, numerically simu¬ 
lated using the finite-element method m- For clarity, 
only the bands with even z—symmetry (out-of-plane di¬ 
rection) are shown, with the green curve corresponding 
to the breathing mode waveguide band. For a lattice 
constant of ao = 480 nm this band has a bandwidth of 
approximately 1 GHz surrounding the 6 GHz frequency 
of the localized acoustic mode of the nanobeam cavity. 
Scaling of the planar waveguide dimensions by scale fac¬ 
tor s' = a/ao, where a is the new lattice constant, re- 
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suits in an approximate scaling of the waveguide bands 
by This is a good approximation for the lower ly¬ 

ing bands of thin structures in which d/ao < 1. As such, 
here we plot the bandstructure in terms of a (unitless) 
normalized frequency s = (//6 GHz)5', where / is the 
physical waveguide band frequency. With the phonon 
waveguide terminated on one end by the nanobeam cav¬ 
ity and on the other end by cutting out a slot in the center 
beam, the continuous waveguide band is transformed into 
discrete waveguide resonances with free spectral range 
/fsr = where Vg = \duj/dk\ is the group ve¬ 

locity of the band and I is the length of the waveguide. 
For those waveguide modes spectrally located within a 
bandwidth je of the localized cavity mode, optical exci¬ 
tation of the optomechanical cavity can be used to detect 
phonons propagating in the waveguide section of the de¬ 
vice. 

Optical excitation and detection of the waveguide 
phonon resonances are performed using a tunable nar¬ 
rowband laser with frequency (ujp) blue-detuned from the 
optical cavity frequency (ujc) by the mechanical cavity 
frequency (cj^), A = oOp — ooc ^ Efficient Stokes 
scattering of the laser pump field occurs for the mechan¬ 
ical breathing mode resonances of the cavity-waveguide 
system which lie within an optical cavity linewidth {k) of 
the bare mechanical frequency For the blue-detuned 
sideband pumping employed throughout this work, the 
mechanical modes also experience amplification by the 
optical pump wave corresponding to a negative (phonon) 
decay rate of 70M = where go is the vacuum 

optomechanical coupling rate characterizing the strength 
of the radiation pressure interaction [2] , Uc is intra-cavity 
photon number determined by the optical pump power, 
and K. is optical cavity damping rate. 

All experiments are performed at room temperature 
and atmospheric pressure. The reflected optical signal 
from the cavity is detected on a photodetector, with 
the beating of the pump wave and resulting Stokes 
waves yielding a photocurrent spectrum containing the 
(mostly) thermal motion of the mechanical resonances of 
the cavity-waveguide system. Figure [TJi shows a typ¬ 
ical microwave photocurrent spectrum of a fabricated 
device with a waveguide length I = 145 /im, cavity-to- 
waveguide coupling rate 7e/27r = 2.4 MHz, and loaded 
optical cavity linewidth K/27r = 0.8 GHz. A series of hy¬ 
bridized modes with free-spectral range /fsr = 4 MHz 
are observed around the main cavity resonance peak at 
co’^/27r = 6.017 GHz, from which one can infer a waveg¬ 
uide group velocity of Vg = 1160 m/s. In order to trace 
out the group velocity evolution along the entire breath¬ 
ing mode waveguide band a set of devices is fabricated 
with the planar dimensions of the waveguide scaled by 
different scaling factors s' = a/ao. The nanobeam cav¬ 
ity parameters are held fixed. In this way the waveguide 
band can be shifted in frequency relative to the nomi¬ 
nally fixed cavity, allowing us to effectively scan across 
the whole waveguide band using the cavity as a frequency 
(and spatial mode) filter. Figure shows the measured 


group velocity evolution along the waveguide band as in¬ 
ferred from the measured /fsr of the waveguide reso¬ 
nances versus normalized frequency s. The green curve 
corresponds to the theoretically calculated group velocity 
curve. We note that in the lower frequency portion of the 
band there exists two /c-points for each frequency, result¬ 
ing in two group velocity curves, a feature also seen in the 
measured data. In analogy to the slow light in photonic 
crystals m, here we measure slow phonon group veloc¬ 
ities [35] near both band edges as well as at the band 
minimum occurring near k ss 0.5(7r/a). The minimum 
measured phonon group velocity is as small as 150 m/s, 
a reduction by a factor of 40 compared to the group ve¬ 
locity of transverse acoustic waves in bulk silicon at room 
temperature. 


III. OPTOMECHANICAL MICROWAVE DELAY 
LINE AND FILTER 

We next study the utility of cavity-waveguide optome¬ 
chanical systems for coherent microwave-over-optical sig¬ 
nal processing. Similar signal processing based on stimu¬ 
lated Brillouin scattering has recently been investigated 
in photonic nanoscale waveguides [384140] . Here we em¬ 
ploy cavities as localized sites for enhancing the interac¬ 
tion of photons and phonons, and phonon waveguides for 
the routing and filtering of microwave signals. The ba¬ 
sic circuit, shown schematically in Fig.[^, consists of two 
optomechanical cavities, labeled L and R, which are con¬ 
nected by a common phonon waveguide. The microwave 
signal encoded in the single sideband of an optical car¬ 
rier with wavelength , is fed into cavity L which coher¬ 
ently drives the localized cavity breathing mode through 
the radiation-pressure force. The excited phonons then 
couple into the phonon waveguide and propagate to cav¬ 
ity R, providing the desired delay and filter functions for 
the microwave signal. In cavity R, the arriving phonons 
modulate the internal optical field, re-encoding the mi¬ 
crowave signal on an outgoing carrier at wavelength Xr. 

A schematic of the experimental set-up used to mea¬ 
sure the two-cavity system is shown in Fig. [^. SEM 
images of the cavity and waveguide elements as they 
are arranged in the optomechanical circuit are shown in 
Fig. An adiabatic fiber-to-chip coupling waveguide 
followed by a ’Y’-branch splitter is used to couple laser 
light into and out of each of the cavity devices using a 
single optical fiber taper. The two nanobeams are de¬ 
signed to have similar mechanical cavity frequencies but 
different optical resonant frequencies so that they may 
be simultaneously and independently probed using two 
separate laser sources (see App. 0 . 

A typical mechanical spectrum of the double-cavity 
system with phonon waveguide length I = 43 gm is 
shown in Fig. H- For this device the two funda¬ 
mental optical cavity modes Ol(^r) have wavelength 
^L(R) = 1539(1559) nm and linewidth = 27r x 

1.15(0.80) GHz. The blue curve corresponds to blue- 


4 



Carrier (A^) 


Carrier (A^) 


Microwave Microwave 
Traveling phonon 


Reciever 


Device 


FIG. 2. Microwave signal processing using optomechanical cavity-waveguide system, a, An optomechanical scheme 
for microwave signal processing relying on phonon routing, b, Schematic of the experimental set-up. VGA, variable optical 
attenuator; FPC, hber polarizer controller; PG, pulse generator; EOM, elector-optic modulator; A-(D)mux, wavelength (de- 
)multiplexer; PD, photodetector; EDFA, erbium doped fiber amplifier; SA, spectrum analyzer; NA, network analyzer; BPE, 
bandpass filter; OS, oscilloscope, c, SEM image of the double-cavity device, showing the ’Y’-branch optical splitter which feeds 
both left and right optical cavities and the intermediary mechanical waveguide that couples the left and right cavities together. 
The mechanical spectrum of the device as measured from cavity L (blue curve) and R (orange curve) is shown in the top inset. 


detuned sideband pumping of cavity L (AL/ 27 r = 6 GHz; 
PpL = 180 /iW; ricL = 1200) and the laser resonant 
with cavity R turned off. The orange curve is for similar 
pumping of cavity R {Aji/27r = 6 GHz; Ppji = 50 /iW; 

= 360). The two spectra are seen to share a number 
of common resonances, a result of the connecting phonon 
waveguide. We label the mechanical resonances as Lj 
(i?/c), with the index j (k) indicating the relative optical 
coupling strength of the resonance to cavity L (R). The 
most intense mechanical resonance peaks in the cavity 
spectra, Li (= R 4 ) and are found to have vacuum 
optomechanical coupling rates of goL^l27r = 0.85 MHz 
to optical mode Ol and ^oi?i/ 27 r = 1.39 MHz to op¬ 
tical mode Or, respectively, determined from the vari¬ 
ation of the mechanical linewidth with pump power. 
The intrinsic Q-factor of the measured mechanical res¬ 
onances are all approximately 1500 as determined from 
the resonance peak linewidths at low optical probe power 
(7i/27r 3.7 MHz). 

We next study the propagation of short phonon pulses 
within the double-cavity system. In this experiment 
the two cavities are simultaneously being probed with 
two separate continuous wave lasers {PpL = 80 /iW, 


ricL = 630; PpR = 50 /iW, Ucr = 360), each tuned to 
the blue sideband of their respective cavity resonance. 
A 20 ns microwave pulse is used to phase modulate the 
laser light input to cavity L at a frequency of 5.952 GHz 
resonant with the Li mechanical resonance, exciting a 
coherent phonon population in cavity L (see App. 0 . 
Figure shows the optically-transduced signal in both 
cavities after the initial excitation pulse. For pulses with 
bandwidth larger than the free-spectral range (FSR) of 
the waveguide modes connecting the two cavities, the 
delay associated with phonons traversing the waveguide 
can be resolved and is given by, r = 1/2/fsr- Multiple 
bounces of the initial phonon pulse are seen in the sig¬ 
nal from each cavity, with a period between bounces of 
58 ns (corresponding to r = 29 ns) consistent with the 
measured FSR of 17 MHz. The pulse propagation can 
be simulated with coupled mode equations (see App. 0, 
and are shown as dashed lines in Fig. |3f- From a fit 
of the model to the measured pulsed response, the ex¬ 
trinsic mechanical cavity-to-waveguide coupling rate is 
estimated to be 7eL(i?)/27r = 35(26) MHz for the L (R) 
cavity, corresponding to an internal efficiency of phonon 
transmission from cavity L to R of 89% (67%) excluding 
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FIG. 3. Pulsed and CW signal propagation, a, A short 
phonon pulse excited by a 20 ns microwave pulse of frequency 
5.952 GHz in cavity L is detected in the two optomechan¬ 
ical cavities as it bounces between them. Dashed lines are 
simulated results, b, GW microwave propagation from cavity 
R to cavity L is characterized by A^-matrix component Srl 
as measured by a network analyzer (see App.[^. Blue curve 
is \Srl\‘^ (in dB) and green curve is the phase of Srl- c, 
Microwave delay at peak of Ri resonance (5.991 GHz) ver¬ 
sus pump power in cavity R (PpR) with fixed pump power 
in cavity L {PpL = 0.2 mW) and microwave signal power 
Ps = —20 dBm. d, \Srl\^ spectrum around Ri resonance cor¬ 
responding to peak delay of 13.3 fis. Red curve is a Lorentzian 
fit. 


(including) phonon waveguide decay. 

Continuous-wave (CW) measurements of the same de¬ 
vice are presented in Fig. in which the frequency (uj) 
of the phase modulation applied to the laser light sent 
into cavity L is varied. The displayed S'—matrix com¬ 
ponent Srl (see App. is for an average optical pump 
power of PpL{R) = 80 /iW for both cavities. Signal de¬ 
lay in this case can be inferred from the phase 0 of Srl, 
r = — where r now includes both the delay asso¬ 
ciated with propagation through the phonon waveguide 
and the resonant storage time in each cavity. Unlike pre¬ 
viously studied optomechanically-induced transparency 
in a single cavity nn, where delay is associated with 


reduced transmission because of the destructive interfer¬ 
ence between direct optical transmission and phonon me¬ 
diated optical transmission, here the delay purely relies 
on phonons as no light directly propagates through the 
waveguide and one can simultaneously realize high trans¬ 
mission efficiency and large delay. For blue-detuned side¬ 
band pumping, which results in parametric amplification 
of the phonon signal, the peak optical transmission gain 
between cavities for the phonon-waveguide mediated pro¬ 
cess can be greater than unity, and is given for excitation 
of mechanical resonance Lj {= Rf^) by. 


G 


j,k,max 


^CL,Cn, 


(1 - Cz,, - CfiJ 


( 1 ) 


where = |7oM,Lj(i?ol/7i is the cooperativity be¬ 

tween the Lj (= mechanical resonance and the OL(i?) 
optical cavity mode. This is the internal gain, and ne¬ 
glects optical loss outside of the cavities and photodetec¬ 
tor inefficiency. 

Figure shows the measured delay for a 5.991 GHz 
microwave signal resonant with the Ri mode as a func¬ 
tion of the optical pump power sent into cavity R for 
fixed pump power PpL = 0.2 mW in cavity L. The delay 
is tunable with optical pump power, reaching a value as 
large as r = 13.3 //s for PpR = 0.21 mW. Figure]^ dis¬ 
plays the measured S'—matrix amplitude spectrum cor¬ 
responding to the longest measured delay, showing a nar¬ 
row 3-dB bandwidth of 17 kHz and a peak internal gain 
of 3 dB. Viewed as a tunable bandwidth microwave fil¬ 
ter [an ED ns, this device realizes an out of band rejec¬ 
tion of up to 70 dB and a thermal-noise-limited sensitiv¬ 
ity of —30 dBm referred to the microwave input of the 
electro-optic phase modulator (see App. 0. 


IV. WAVEGUIDE-MEDIATED DISTANT 
MECHANICAL COUPLING 


In addition to guiding the propagation of acoustic 
waves, phonon waveguides may also be used to medi¬ 
ate coupling between distant optomechanical cavities via 
virtual phonons, i.e. without energy distribution in the 
mediating waveguide modes. Such a scenario can be 
used to realize highly efficient signal conversion between 
cavity elements, or in the case of mechanical sensing, 
may be employed to realize new sensing modalities in 
which read-out and sensor are separated spatially. As 
recently demonstrated for optical photons in photonic 
crystals m, when the cavity-to-waveguide coupling rate 
is substantially smaller than the FSR of the waveguide 
modes, 7 e/ 27 r ^ / fsr , the excitation of real waveg¬ 
uide phonons is suppressed and only second-order virtual 
transitions between cavity and waveguide modes con¬ 
tribute to the inter-cavity coupling. For nearly degen¬ 
erate bare cavity modes the waveguide mediated inter¬ 
cavity coupling strength scales as % = ■sJ'^eL^eR, and for 
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FIG. 4. Waveguide-mediated distant mechanical cou¬ 
pling. a Mechanical spectrum of a device with two distantly- 
coupled cavities as measured from cavity L (blue) and R (or¬ 
ange). Modes and C- are the (approximately) even and 
odd hybridized mechanical cavity modes with the simulated 
modal profile shown in the inset. The weakly transduced Wi 
and W 2 resonances are waveguide modes, b Measured phonon 
Rabi oscillation between the two distant cavities. 


7e > 7i+7oM such a configuration can provide both high 
internal (i.e. phonon routing between cavities) and exter¬ 
nal (i.e. phonon-photon exchange in the cavity) efficiency 
for the microwave filter demonstrated in the previous sec¬ 
tion. 

Figure Ef shows the measured (thermal) mechanical 
spectra of a device with two optomechanical cavities cou¬ 
pled together by a 30 /im long phonon waveguide which 
is designed to have 7e/2/pi ^ /fsr- This device is simi¬ 
lar to that in Fig. with one key difference: the central 
holes of the waveguide have been removed, reducing the 
cavity-to-waveguide coupling. An optically-transduced 
spectrum measured through both the left (L) and right 
{R) optical cavities is shown. The two hybridized me¬ 
chanical cavity modes, labeled (7+ and (7_, have a fre¬ 
quency splitting _/27r = 4 MHz and are separated 

from the nearest waveguide mode by (5 = 17 MHz. Both 
of the (7+ and (7_ mechanical resonances are measured 
to have large optomechanical couplings of over 700 kHz 
with both optical cavity modes, indicating strong mixing 
of the L and R mechanical cavity modes. The weakly 
transduced waveguide modes are measured to have a FSR 
of /fsr = 54 MHz. 


From the matrix of optomechanical rates between the 
hybridized mechanical modes and the (uncoupled) op¬ 
tical cavity modes we determine a waveguide-mediated 
coupling of H/27r = 1.94 ±0.06 MHz and a bare mechan¬ 
ical mode splitting of /S.lrI2t: = 0.98 ± 0.48 MHz (see 
App. |§. The coupling rate between nearly degenerate 
L and R mechanical cavity modes can be shown to be 
given to lowest order by H = 7e/ sin(7r^//FSR), yield¬ 
ing an average cavity-to-waveguide coupling of 
1.62 MHz, consistent with the numerically-simulated de¬ 
sign value for this device. Despite the distant separation 
between the two cavities, phonons can nevertheless tun¬ 
nel between them via virtual excitation of the waveguide 
modes. To show this, a 50 ns pulse of a 5.95 GHz mi¬ 
crowave modulation signal is applied to the optical laser 
beam input to optical cavity L, locally exciting a coherent 
phonon population in the left cavity. Following the pulse, 
Fig.[^ shows the resulting oscillation of the phonon pop¬ 
ulation between the two cavities. The measured period 
of oscillation is 250 ns, corresponding to a Rabi frequency 
of 4 MHz commensurate with the hybridized mode split¬ 
ting _ /27r. This result should be contrasted with the 

previously demonstrated phonon pulse bouncing shown 
in Fig. where the period is set by the free spectral 
range of the waveguide modes. 


V. DISCUSSION 

The devices presented in this work, involving pairs of 
optomechanical cavity elements connected together by a 
phononic waveguide, demonstrate the feasibility of creat¬ 
ing small optomechanical crystal circuits for microwave- 
over-optical signal processing and novel sensing modal¬ 
ities. Advantages of the optomechanical approach are 
myriad. For instance, the use of radiation pressure as op¬ 
posed to the piezo-electric effect to couple to microwave 
phonons allows one to work with a wider variety of ma¬ 
terials, such as silicon in the devices studied here, and to 
realize high efficiency of signal conversion. Purely elec¬ 
tromechanical nanoscale devices [43l[44], by comparison, 
suffer from capacitance and impedance mismatch which 
limits efficiency. Implementation of two-dimensional op¬ 
tomechanical crystals, as has been recently demonstrated 
in Ref. |45] , can also be used to create multi-tap phonon 
waveguide structures capable of realizing diverse spectral 
filter shapes and dispersion control without the need for 
additional electrical or photonic hardware. 

Further scaling of optomechanical circuits, to arrays 
of a large number of cavity elements, is also interest¬ 
ing. A straightforward one-dimensional extension of the 
currently studied structures to a chain of mechanically- 
coupled optomechanical crystal cavities, for instance, can 
be used to study mechanical collective dynamics m- 
An even richer set of physical phenomena and device 
functionality, however, may be realized in optomechan¬ 
ical circuits in which both phononic and photonic con¬ 
nectivity of cavity elements are employed. For exam- 
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pie, building on the tight-binding mechanical coupling 
of cavities demonstrated in this work, an acoustic and 
optical isolator may be realized using a phase-correlated 
optical pumping scheme in which the cavities are also 
coupled together by a photonic waveguide [IQIST]. In 
this scheme [48], analyzed in App. the interference 
between optical and mechanical hopping amplitudes be¬ 
tween the cavities picks up the phase difference between 
the laser pump beams at each cavity site, giving rise to 
the required non-reciprocal behavior. With the help of 
post-processing methods, such as scanning probe nano¬ 
oxidation lithography m which has been successfully 
employed to tune photonic crystal cavities m, the task 
of aligning both mechanical and optical resonances of the 
optomechanical cavities could be fulfilled. This would al¬ 
low, at least in moderate lattice sizes, the study of topo¬ 


logical phases and effective magnetic fields for photons 
and phonons [3QH3^ . 
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Appendix A: Methods 


3. Pulse and CW measurements 


1. Fabrication 

The devices are fabricated using a silicon-on-insulator 
wafer with silicon device layer thickness of 220 nm and 
buried-oxide layer thickness of 2 /im. The device geom¬ 
etry is defined by electron beam lithography followed by 
inductively-coupled-plasma reactive ion etching to trans¬ 
fer the pattern through the 220 -nm silicon device layer. 
The devices are then undercut using an HF:H 20 solu¬ 
tion to remove the buried oxide layer, and cleaned using 
a piranha etch. 


2. Experimental Set-Up 

Two lasers are used to probe the two nanobeam 
cavities individually, one of which can be encoded with 
microwave signals by an electro-optic phase modula¬ 
tor. The output of the two lasers are combined in a 
wavelength multiplexer and coupled into the device 
through a tapered optical fiber. The on-chip optical 
waveguide is split into two paths to transmit light 
into the two nanobeam cavities through side-coupled 
optical couplers. The reflected light from the cavities 
is re-collected by the optical fiber and passed through 
an optical circulator to separate it from the input light. 
A fiber-based wavelength de-mult ip lexer with 20 nm 
channel spacing is then used to separate the two cavity 
signals. Finally, the reflected light is amplified by an 
erbium-doped fiber amplifier and sent to a high speed 
photodetector (12 GHz bandwidth). The electrical 
signal from the photodetector is measured either with 
a spectrum analyzer, network analyzer, or high speed 
oscilloscope. 


For the phonon pulse measurement, a pulse genera¬ 
tor is used to gate the microwave source which drives 
the electro-optic modulator with a repetition rate of 
20 kHz. The photodetected optical signal is mixed down 
to 300 MHz with a second microwave source which is 
locked to the one driving the electro-optic modulator. 
The down-converted signal is filtered, amplified and av¬ 
eraged 10 ^ times on an oscilloscope with sampling rate of 
10 GSa/s. For the GW measurements a network analyzer 
is employed, with the microwave output of the analyzer 
(port 1 ) connected to the electro-optic modulator used 
to phase modulate the light sent into cavity L, and the 
photodetected optical signal reflected from the cavity R 
is connected to port 2 of the analyzer. We define the S- 
matrix component Srl = § 21 , where .§21 has been nor¬ 
malized to remove the effects of optical insertion loss into 
the devices, optical pre-amplification, and photodetector 
gain. 

Appendix B: Simulation of phonon pulse 
propagation 

In this section, we show propagation and bouncing of 
phonon pulses in the cavity-waveguide system (Fig. 3a) 
can be well simulated by a group of coupled mode equa¬ 
tions using input-output formalism. The dynamics cap¬ 
tured by the coupled mode equations is a phonon pulse 
travelling in a waveguide terminated by two cavities with 
bare mechanical frequency uJmL,R and waveguide cou¬ 
pling rate We approximate uJmL,R to be the fre¬ 

quency of cavity-dominated modes Li and Ri in the sim¬ 
ulation. Since the response time of the optical cavity is 
much shorter than that of the mechanical cavity, we can 
exclude the dynamics of optical modes from these equa¬ 
tions. Thus, the coupled mode equations can be written 
as follows. 


^ - i<7oLaoLa;Le-“**0(t)e(r - t) + (Bl) 

^ (B2) 

bin,L{t) (B3) 

binMt) = e~°‘\^/^bL{t - tw) - bin,L{t - tn,)), (B4) 


where a^L and are the amplitudes of optical pump 
and its red sideband in the left cavity, r is the duration of 
excitation pulse, ujs is the frequency of pulse, 0 (t) is the 
Heaviside step function, 7 is the effective decay rate of 
the excited mechanical mode, a ^ j/vg is the waveguide 
loss rate, and = 1/(2/fsr) - l/(7eL +7) “ l/(7eii+ 7) 
is the single trip time the pulse spent in the waveguide. 


From the mechanical spectrum we And 7 = 27r x 2.1 
MHz for Li mode (the main coherently-driven mode) 
during the pulse measurement; and by fitting the pulse 
tails detected in each cavity we And 7 eL = 27r x 34.7 MHz 
and 7 ei? = 27r x 25.5 MHz. Using these parameters, |5 l| 
and \bR\ can be numerically calculated from the coupled 
mode equations and the proportional voltage signals are 
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TABLE L Mechanical mode parameters 



5o/27r (MHz) 

jil 271 (MHz) 

(MHz) 

Li 

0.85 

3.7 

34.7 

L2 

0.75 

3.7 

L 3 

0.63 

3.9 

Ri 

1.39 

3.6 

25.5 


shown in Fig. 3a. The simulated result captures the main 
features of the measured pulse data. In particular, the 
pulse splitting observed from cavity R is due to the fact 
that the pulse frequency is not in resonant with cavity R 
and thus experiences destructive interference inside this 
cavity. 

The phonon transfer efficiency from cavity L to cavity 
R is about e“'T/(2/FSR) ~ 67%. The phonon transfer ef¬ 
ficiency from cavity to waveguide for cavity L and R is 


7eL(i?)/(7eL(i?) +7) ~ 94%(92%). 

We summarize the measured mechanical mode param¬ 
eters of the device in Fig. 3 in Table [Tj where go for the 
Lj{Rk) modes is with respective to OL(i?) optical cavity 
modes. 


Appendix C: Derivation of the microwave S'—matrix 
for the optomechanical cavity-waveguide system 

Here, we derive the S'—matrix for a microwave signal 
traversing the optomechanical cavity-waveguide system. 
We assume the mechanical amplitude is small such that 
only the first-order optical sideband needs to be consid¬ 
ered. In the next section, we will verify the small me¬ 
chanical amplitude assumption. 

The Hamiltonian of the system under continuous wave 
operation involves two optical cavity modes aL^R with 
frequency 0 JcL,r parametrically coupled to a common me¬ 
chanical mode 6, 


E huJck^\ak + huJrnb^b + ^ + E ihy/K^apke - &{), (Cl) 

k=L,R k=L,R k=L,R 


where the last term is the pumps of the two optical cav¬ 
ities. For simplicity we assume the pumping lasers are 
blue-detuned from the cavity resonances which is true 
for all of our experiments. Suppose the pumping laser 
for cavity L is modulated at frequency cj by a microwave 
signal, then the operators of the system can be decom¬ 
posed into carriers and sidebands, 

ak = b = (C2) 


where we only keep the red sideband of the pumping 
lasers because of rotating wave approximation, given the 
sideband resolved condition ujm ^ kik of our device. Sup¬ 
pose the pumps are strong enough such that the car¬ 
rier operators can be treated as static variables, then the 
equations of motion of the system can be derived after 
substituting Eq. |C2| into Eq. EH 

'i'gOkk^OkP— \/ ^e/c0^in,k: 

(C3) 

-iup- = -{iuJm + “ E *5ofeaofeatfe’ (^4) 

k 


where Afc = copk-uick ~ 
in the frequency range \uj — uj 


Solving Eq. C3 and Eq. C4 
ml ^ we obtain 


P- = 






-a 


in,L5 


(C5) 


— —\/KeROt+R 

_ 4:goLgoR^KeLReR/{KLKR)aoi^aoR 

From Eq. |C6[ peak optical gain at w = Um is 


(C6) 

0^in,L • 


_ logout,it:P _ ^ClCr 
|<^in,Lp (1 - C'l - Cr)‘^ ’ 


(C7) 


where = \lOM,L{R)\/li is the coorperativity of me¬ 

chanical mode b with optical modes 

Using the result of Eq. |C6[ the microwave signal trans¬ 
fer S'—matrix can be derived 


Srl 

Una,] 


(C8) 


Una, out 

^oL^oit:GeGEDFA(^^cit:^m/ l^eR)o^o\it,R(y^R 
(2f4/7r) ((ain,L/ 

^oL^oi?GeGEDFA ^gOL 90 RblOcRio‘^/ i?) |ctOL P 




2WA 






‘^9lkWok\‘^ 


where Una, out and Una, in are the output and detected 
electrical voltage of the network analyzer respectively, 
VoL,R is the optical loss of the input and output ports of 
the device and fiber respectively, Gedfa and Ge are the 
gain coefficients of EDFA and photodetector respectively, 
and Utt is the voltage required to produce a phase shift 
of TT of the electro-optic modulator. 
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Appendix D: Analysis of linear operation and noise 
characteristics of the optomechanical microwave 
filter/delay line 

In this section, we examine the assumption of weak 
mechanical amplitude under strong optical pump and an¬ 
alyze the optomechanical microwave filter/delay line per¬ 
formance in terms of linearity and noise characteristics. 
We find that thermo-optic effect constrains the mechan¬ 
ical amplitude due to saturation of the optomechanical 
gain. This effect sets the linear operation range and the 
suppression of the mechanical thermal noise. 

The thermo-optic effect induced optical cavity fre¬ 
quency shift can be described by the following equa¬ 
tions [51] 


where ric is cavity photon number, nsi, r, and is the 
refractive index, thermal resistance, and two-photon ab¬ 
sorption coefficient of silicon respectively, c is the speed 
of light, Vtpa is the cavity volume for two-photon ab¬ 
sorption, and A is a perturbation t heor y coeffic ient A = 


|E(r)pdr 


:. Substituting Eq. 


D2 


into Eq. 


/nsi(To)2|E(r)|2dr- ___ 

using the parameters of silicon given in Ref. along 
with A ^ 7.5 X 10“^, Vtpa ~ (A/nsi(To))^, we have 


D1 


and 


6ujc = <^^05 ^ ^ —33.9 Hz. 


(D3) 


6ujc = -cOcnsiiTo) 


5T = 


dnsijTo) 

dT 


r^c 




AST, 

(Dl) 

nl, 

(D2) 


TPA 


W e pr oceed to include the term of therm o-op tic effec t 
(Eq. |D3[ ) to the equations of motion (Eqs. |C3| and C4), 
which are then modified to be 


iu:a+ = (z(A + ^(|q;o|^ + |a+|^)^) - -Wo^oldl, 
-iujp- = -{iuim + - igoOiQa\ - An = 


(D4) 

(D5) 


where we have denoted and P- as the static value of 
the corresponding operators without input optical side¬ 
band signal and we have included the mechanical thermal 
noise input. Also, we specifically consider the operation 
with the Ri mode, and thus ignore optical cavity L which 
has much weaker coupling with Ri compared to optical 
cavity R. Eqs. |D4| and |D5| can only be solved numerically 
for a generic pump condition. To reveal the thermo-optic 
effect on the mechanical amplitude, we consider a special 
pump condition corresponding to the original threshold 
of mechanical self-oscillation, i.e. |aoP/(/^7i) = 1. 
In this case, we can analytically solve for the down- 
converted photon number and mechanical amplitude at 
uj = ujjn from Eqs. |D4| and D5 , assuming A = and 
|d+| > l^ol, 


I- 12 _ 


l«+P = (- 


4^2 


7i 44 


^th • 


(D6) 

(D7) 


Eor A = 27r X 0.8 GHz, 7 ^ = 27r x 3.6 MHz, qq = 2tt x 1.39 
MHz, nth = ~ 1000, we have Id+p 1900, 

1^0P ~ 380, and \P-\^ ~ 3.7 x 10^. It is the optical 
resonance shift induced by thermo-optic effect that sat¬ 
urates optomechanical gain and prevents runaway of the 
mechanical amplitude at the threshold. 


Now we can estimate whether the mechanical ampli¬ 
tude is large enough to induce nonlinearity through exci¬ 
tation of higher order optical sidebands. The nonlinear¬ 


ity arises due to pump saturation and occurs when the 
amplitude of the first order optical sideband significantly 
deviates from being linearly proportional to the mechan¬ 
ical amplitude, i.e. approximation Ji{z) ^ | breaks 

down [ 52 ], where 2 ; = + ‘^/^m is the normal¬ 

ized mechanical amplitude. Eor |/d_p ~ 3.7 x 10^ calcu¬ 
lated above, the deviation is only about 1%. In the ex¬ 
periment, we find for the largest pump power PpL,R ~ 0.2 
mW, |/d_p « 2.0 X 10^, which gives 2 ; = 0.18 and a lin¬ 
ear deviation of 0.4%. As a result, scattering into higher 
order optical sidebands does not need to be included in 
Eq. (|D4|) and ([D^. 


Next, we consider the response of the mechanical os¬ 
cillator to a small input optical sideband signal by per¬ 
turbative expansion of Eqs. |D4| and |D5l In this case the 
coherent mechanical amplitude of Eq. |C5| is modified to 
be 


P- = 


'ig0L^/l^eL-^O^0L 


W 2 l^k -2i5 


^in,L5 


(D8) 




where 5k = ^(|(ao/cP + is the thermo-optic-effect 

induced optical frequency shift. Eqs. |C6| an d |C8| can be 
modified correspondingly. According to Eq. |D8[ the ef¬ 
fective mechanical loss rate is 7eff = 'll — ^k ^ 

^9ok\^ok\ rppg deviation from a linear re- 

spouse can be caused by the additional cavity photons 
from the input signal, and is characterized by the ra- 
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Frequency offset (kHz) 

FIG. 5. Phase noise suppression ratio between microwave 
signal power of -15 dBm and -30 dBm at largest optical pump 
level. Solid line is theoretical value calculated using Eq. |D11| 
and dots are experimental data. 


tio r = |o+/cP/|q;+/cP (in our device the contribution 
is mainly from At the theoretical self-oscillation 

threshold, we find the 1 dB compression point of the 
S'—matrix to be equal to a microwave power of -19 dBm 
(assuming |(aoLp = experiment, we find 

for the largest optical pump power PpL,R = 0.2 mW 
(which is slightly above the self-oscillation threshold), 
the 1 dB compression point occurs at a microwave signal 
power of -15 dBm. For reversed operation (cavity R as 
input), the 1 dB compression point is reduced by a factor 
of 70M,L/70M,ii (assuming |q;ol|^ = |aoiiP)- 

We now analyze the noise characteristics of the op¬ 
tomechanical cavity-waveguide system. The dominant 
form of noise is from thermally excited phonons in the 
system. From the thermally-added mechanical noise ref- 
ered to the input signal (/^e/^)~^7i^th/|7oM| [23], we 
define the noise-equivalent optical signal power 

Pne = —fiwc—27rP, (D9) 

70M 

where B is the bandwidth of the coherent signal and all 


the quantities are referred to the input cavity. Then the 
noise-equivalent microwave signal power is 


T/2 _ 

^NE — 


4 -Pne .^2 

7r2 Pp - 


(DIO) 


We find that for the largest pump power when operating 
at Ri resonance, if the input port is cavity L, the noise 
equivalent microwave power is -30 dBm; if the input port 
is cavity i?, the noise equivalent microwave power reduces 
to -70 dBm because of the significantly enhanced 70 M of 
cavity R with Ri mode. 

For a self-oscillating mechanical oscillator, the intrinsic 
oscillator noise can be suppressed by the injection of an 
external coherent signal [53] . The suppressed phase noise 
(ignoring input signal noise) can be modeled by [54] 


LUm) — ^ iOm), (Dll) 


where is the intrinsic phase noise spectral density 

without injection and p = is the ratio between 

injected mechanical amplitude and free-running ampli¬ 
tude. Experimentally, we infer the phase noise from the 
measured noise power spectral density using the defini¬ 
tion ^'^^(ci;) = Si)b(u;) — f Sbb{^)duj. At the largest optical 
pump level, the phase noise suppression ratio between the 
microwave signal power -15 dBm (1 dB S'—matrix com¬ 
pression point) and -30 dBm (noise equivalent power) is 
shown in Fig. The model of Eq. |Dll| explains well the 
measured noise suppression level in the offset frequency 
{uj — ujrn) range between 7eff/4 and 37eff/4 (yeff = 27r x 17 
kHz). 


Appendix E: Waveguide-mediated cavity-cavity 
coupling 

1. Analytical derivation 

When 7 eL,i?/ 27 r ^ /fsr? and all the waveguide modes 
are outside of the line width of the cavities, two mechan¬ 
ical cavities can acquire a waveguide mediated coupling. 
The generic Hamiltonian describing this case is 


H = hWmLh^J)L + hujkb\hk + KgLkb^bk + gtkbLb\) + T] ^^dRkb^Rbk + 9*RkbRb\), (El) 


where b\{bk) is the creation(annihilation) operator of the 
k-th waveguide mode, gLk = V^tEJfsr and guk = 
(—)^\/ 27 ei?/FSR are the coupling coefficients of the left 
and right mechanical cacity modes with the k—th waveg¬ 
uide mode, and the summation is over all waveguide 
modes. Note the relative sign between gLk and gRk comes 
from the symmetry of waveguide modes with respect to 


the center of the two cavities. 


In the degenerate cavity case, i.e. we 

can calculate the coupling V between the two cavities 
mediated by the waveguide modes using second order 
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perturbation theory as follows, 


27r 


E 


9Lk 9Rk 
277 277 

^— ^/fsr 


^/leLjeR 
27r 


E 


i-Y 


k= — oo /fsb 
^/leLleR 1 
27r sin-y— tt’ 

JFSR 


kn 


(E2) 


where (5 is the frequency difference between the degener¬ 
ate cavity modes and their nearest waveguide mode. For 
the perturbation theory to be valid, we require 


\ 9 L,Rk |/27r = \/27eL,R/FSR/27r < 5. (E3) 


We consider two special cases. First, when S = /fsr/2, 
i.e. the two cavity modes are in the center of two nearest 
waveguide modes, then V = In the other limit 

when S < /fsr, then V ^ ^y^eL7eR^^• One can prove 
in this case, V < ^/2JeL,RfFSR = \gL,Rk\^ by taking into 
account of the condition of Eq. |E3[ Different from the 
previous case, here the contribution to the coupling can 
be almost exclusively attributed the nearest waveguide 
mode (/c = 0). 


2. Measurement modeling 

We model how to estimate the phonon waveguide me¬ 
diated coupling between two optomechanical cavities. In 
the device for demonstrating waveguide-mediated cavity 
coupling, we used a phonon waveguide without air holes. 
The band structure of the waveguide is shown in Eig. [^. 
Such a waveguide provides large phonon group velocity 
and thus large free spectral range in order to isolate cav¬ 
ity modes from waveguide modes. 

In the sample device shown in Eig. 4 of the main text, 
virtual phonons in waveguide mixes the mechanical cav¬ 
ity modes Ml and Mli of cavity L and R into hybridized 
modes and C -. The coupling strength can be inferred 
by measuring the optomechanical coupling of hybridized 
modes and C- with the two optical modes Ol and 
Or based on the following model. (7+ and C- can be ex¬ 
pressed as linear superposition of Mr and Mr^ assuming 
ignorable energy distribution in the waveguide (as 
are well separated from waveguide modes), 

Cr = aiMR + a2MR^ (E4) 

C-=PiMl+P2Mr. (E5) 

The superposition coefficients ai and /3i satisfy the fol¬ 
lowing relation 


Pi I _ 1^1 _ 2|E| 

'/32'“ + A + 


(E6) 


where V is the waveguide mediated coupling between Mr 



FIG. 6. a Band structure of shielded phonon waveguide with¬ 
out air holes. Inset is the modal profile with k — 0.87r/a. 
b,c Characterization of optomechanical coupling between hy¬ 
bridized mechanical modes (7+ and C- with the optical modes 
of cavity L (b) and R {c). 7 = 7 z+ 7 om is effective mechanical 
decay rate and ric is the cavity photon number. 


and Mr^ ^lr is the fabrication induced fr equency dif- 
ference between Mr and Mr, and A+_ = \/ ^\r + 41/^ 
is the frequency difference between the hybridized modes 
Cr and (7_. On the other hand, the ratio of the super¬ 
position coefficients is directly related to the measurable 
optomechanical coupling 


9 +l _ 1^1 g+R _ 1 ^ 
9-R~^p2 
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where, for example, stands for the coupling between 
mechanical mode and optical mode Ol- 

We inferred the optomechanical couplings by mea¬ 
suring the pump dependent effective mechanical damp¬ 
ing rate of (Fig. c) using the relation 7 = 

T'i “ obtained 

g+L = 927 ± 25 kHz, g_L = 710 ± 55 kHz (E 8 ) 
g+R = 683 ± 90 kHz, g_R = 870 ± 35 kHz (E9) 

As a result, the ratio of superposition coefficients is 

1^1 = 1.30 ±0.14, 1^1 =0.79 ±0.14. (ElO) 

Pi P2 

Along with _/27r = 4.0 MHz as read from the spec¬ 

trum, we find, according to Eq. |E 6 | 

|E|/27r = 1.94 ± 0.06 MHz, ALi?/27r = 0.98 ± 0.48 MHz. 

(Ell) 

To compare with the analytical formula of V (Eq. |E 2 [), 


in this device, we have /fsr = 54 MHz and S = 17 
MHz. And using the measured E/27r = 1.94 MHz, we 
find y7eL7eR/27r = 1.62 MHz. 


Appendix F: Optical non-reciprocity based on 
distantly-coupled optomechanical cavities 

We provide theoretical analysis of achieving optical 
non-reciprocity in the distantly-coupled optomechanical 
cavities, and show its viability based on the typical pa¬ 
rameters of our fabricated devices. In this case, the 
waveguide connecting the two cavities should support 
both guided mechanical and optical modes, and the two 
optomechanical cavities are designed to be identical. As 
already been demonstrated in our experiment and in 
Ref. [20] respectively, such a waveguide can mediate a 
tight-binding type coupling for the mechanical and op¬ 
tical cavity modes. The Hamiltonian of the system can 
thus be written as follows. 


E hjJcka\ak ± J{a\aR ± ULaJj) ± E h^mkh\hk ± V{h{hR ± (El) 

k=L,R k=L,R 

± E f^9k{bl + bk)alak + ^ - a\), 

k=L,R k=L,R 


where J and V are the waveguide mediated optical and 
mechanical coupling strength, and the last two terms are 
the optical pumps in the two cavities which have a same 
frequency and correlated phases. 

The optical non-reciprocity arisen from this system 
can be intuitively understood from a schematic shown 
in Fig. The input optical signal undergoes a Mach- 
Zehnder type of interference through the system: one 
path is the direct photon hopping and the other path is 
through radiation-pressure interaction induced transition 
to phonon and phonon hopping. The phase of the lat¬ 


ter path involves the phase difference of the two pumps, 
which is ipL — ^R for one direction and (fR — for the 
reversal direction. Such a non-reciprocal phase resem¬ 
bles an effective magnetic flux for photons, resulting in 
the non-reciprocal transmission [48] . 

We first consider the case when both cavities are be¬ 
ing pumped with blue detuned lasers {ujp — uOck = ^mk)- 
The equations of motion of the system, in terms of red 
optical sidebands of the pumps, can be derived from the 
Hamiltonian of Eq. |E1| using rotating wave approxima¬ 
tion, given sideband resolved condition ujrnk ^ 


daL 

dt 

daR 

dt 


= {'idh -~ iJaR - - y/ReLdL.in, 

(E2) 

= - — iJaR — igROLRd'^^dR - 

(E3) 

= -{iuJmL + ~ '^VdR - igLOihd'^^oTl, 

(E4) 

= -{iuJrnR + y ~ “ WROiRd"^^O^ r, 

(E5) 

1 

where 6 ^ = oOp — uuck and 

is the steady state am- 


plitude of the local optical cavity mode, which is related 
to the pumping amplitudes as follows 
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FIG. 7. a Schematic of the optical non-reciprocity in distantly-coupled optomechanical cavities stemmed from the non-reciprocal 
optical pump phases and interference between photonic and phononic transmission paths, b-e Transmission coefficient b(d) 
and isolation ratio c(e) for blue (red) detuned optical pumps with practical device parameters 7^ = 27r x 3 MHz, ql^r = 27rx0.8 
MHz, cavity photon number ncL,R = 2000, rl,r = 27r x 2 GHz, ReL,R = 27r x 1 GHz, (^l — 4>r = f , J = 27r x 422 MHz, and 
y = 27r X 10 MHz. 


- kr(l)/ 2)^KeL(R)apL{R)e +iJ^/K^^L)apinL)e 

{iSL - K.L/‘^){iSR- k.r/2) + R 


(F6) 


We find the steady state amplitude is approximately 
Under the condition ~ ujmk ^ 
Ji which means each cavity is effectively only driven 
by its own optical pump. This can be intuitively under¬ 
stood by the fact that even and odd hybridized optical 
cavity modes are driven equally (as 5^ ^ J) and thus the 
amplitude of one local cavity mode is cancelled out and 
effectively not being driven. Thus, each cavity-enhanced 
optomechnical coupling can be independently controlled 
by the pump. 


After solving the equations of motion, we calculate the 
ratio between right transmission coefficient Tr and left 
transmission coefficient Tr of the optical signal and find 


Tr 

Tl 


J - 


VgLgROCLOtR 




pi{(t>L—4>R) 


J 


VgLgROtLOtR 


(i (u; -I. ) + ^ ) R) + ^ ) +^ 


^ — i{4>L—4>R) 


(F7) 

Interestingly, this ratio is not explicitly dependent on 5^ 
and Rk us an intrinsic pro perty of the device. At the 
poles w = {u)mL + i^mR ± vTwTFTF+F 2)/2, i.e. 
frequency of the hybridized mechanical modes, and as¬ 


suming t we have 

J + n VgLgROtLOtR ^i((t)T,-cl)n) 

Tr _ _ + _ .pgx 

Tl ~ J^i _ VgLgRotLotR • v ) 

Thus perfect non-reciprocity, i.e. one direction has van¬ 
ishing transmission while the other direction has maximal 
transmission, can be achieved by satisfying the following 
condition 


0 L 


<t>R = ±|, 


_ VgLgRO^LO^R _ 

+ {^mL ^turY 


(F9) 

Under this condition, the transmission coefficient for the 
through direction is (for simplicity assuming ujmL = 

^mR) 


— \/ l^eLl^eR / 


2 gLgROtLOtR 

_ Ji _ 

gLgROtLOtR \( I^R _ gLgROtLOtR \ • 
li A 2 li ^ 

(FIO) 


Similar results can be derived for the case of red de¬ 
tuned pumps, and we find the perfect non-reciprocity 
condition (Eq. |F9[ ) remains the same; while the trans¬ 
mission coefficient for the through direction at poles is 
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given by 

2 QLgRC^LaR 

— ^/l^eLl^eR , gLgRaLO^R \( i^b. i gLgRaLO^R \' 

\ 2 ^ 7i A 2 7, ) 

(Fll) 

We note, in general, an amplified transmission in blue 
detuned case and an attenuated transmission in red de¬ 
tuned case for the through direction at poles. One 



Based on our fabricated devices, real izing the condi¬ 
tions for perfect non-reciprocity (Eq. |F9| ) is quite promis¬ 


ing. For a typical device with 7 ^ = 27r x 3 MHz, Ql = 
qr = 27r X 0.8 MHz, maximal available cavity photon 
number Ucl = '^cR = 2000 , and assuming UmL = ^mRi 
Eq. |F9| determines Jmax = 422 MHz. Thus, as long as 
J < -^max in this device, perfect non-reciprocity condi¬ 
tion can always be achieved by tuning the pump power 
and phase. We have numerically simulated waveguide 
(without acoustic shielding) mediated optical coupling 
between two optomechanical cavities to be 500 MHz and 
less. Thus it is indeed viable to demonstrate optical non¬ 
reciprocity in our devices. Using these parameters along 
with 5k = ojmki = 27r X 2 GHz, R^k = 27r x 1 GHz and 
(t)L - (t)R = J = 2 t: X 422 MHz, U = 27r x 10 MHz, 
we plotted the transmission coefficients and isolation ra¬ 
tio as calculated from Eqs. |F2f[F5| for the blue-detuned 
pump In Fig. [^,c and red-detuned pump In Fig. [ 3 i,e, 
respectively. 













